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We discuss a laser-trapped cold-atom superfluid qubit system. Each qubit is proposed as a 

macroscopic two-state system based on a set of Bose-Einstein condensate (BEC) currents circulating 

in a ring, cut with a Josephson barrier. We review the effective low energy description of a single 

BEC ring. In particular, it is demonstrated that such system has a set of metastable current states 

— ■ which, for certain range of parameters, form an effective two-state system, or a qubit. We show how 

^^ ' this qubit can be initialized and manipulated with currently available laser-trapping techniques. 

, We also discuss mechanisms of coupling several such ring qubits as well as measuring individual 

O ■ qubit-ring systems. 

Q 
o 
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Since the first successful experiments on Bose-Einstein condensation of alkaU gases, optical cooling and trapping 
jy^', have become a standard experimental tool to study quantum degenerate gasses at ultra-low temperatures. The laser 
C^ ' trapping of neutral atoms has been rapidly developing over the past decade. State-of-the-art traps are capable of 
^^iJ] forming precise time- varying potentials with micron-size resolution [l|. High densities (~ 10^'^"^^ cm""^) and wide 
range of interactions have been achieved at nK temperatures [2|, 3 ■ Not surprisingly, cold atom systems have become 
5 ' a perfect tool to test fundamental physics from basic quantum mechanical principles and quantum information to 

complex strongly correlated quantum states of matter (j-Q ■ 

^H' Trapping of neutral atoms is based almost entirely on resonant scattering of light. It utilizes the internal energy 

; , structure of an atom to the control radiation pressure of the trapping (cooling) laser beam. An accurately configured 

d ' laser field creates an effective (attractive or repulsive) potential profile for each atom, see @, [lOhlSj for further details. 

Each atom is indistinguishable and has either integer or half-integer spin depending on its nuclear and electronic 

T ■ content. Hence, at the temperatures of typical cold-atom experiment, atoms obey Bose or Fermi statistics, e.g. 

^ _ typical "bosons" are ^^Rb, ssRb, '^Li, ^^^Cs, etc. 

Quantum information has been introduced in a cold-atom system on various levels |14l . [l5j. It is natural, for 
example, to define a qubit (two state system) via two internal states of an atom. In this case, each atom records a 
single bit of quantum information. This approach, however, requires each atom to be addressed separately. A similar 
problem appears when the qubit is introduced via a set of spatially localized states (e.g. in adjacent wells of an optical 
lattice potential) of an atom or a dilute Bose-Einstein Condensate (BEC). The complication is due to the fact that 
[^ ■ the number of atoms, N ^ in a BEC (cold-atom) experiment fiuctuates significantly, about 10% (or at least as viV), 
QQ ' from run to run. Therefore any qubit system dependent on the number of atoms becomes problematic. 
^^ Another way to introduce a two-state system is via a collective phenomena. In this case the qubit has to be formed 

'^ [ by a pair of distinct macroscopic states that are sufficiently far away from other states of the multiparticle system and, 
\f-\ . at the same time, have the energy difference between these lowest states small enough to allow measurable dynamics 
f^ ' [16|. Already in the early stages, the experiments on BEC of trapped atoms were focused on the dynamics in a 
^D , double- well(dot) trapping potential [ij]- However, despite the appealing similarity with a single-particle two-state 
^"^ ■ system, the BEC system of such geometry presents no easy way to achieve superposition of macroscopic states with 
^ , measurable two-state dynamics. Consider, for example, a typical system of non-interacting (weakly interacting) boson 
atoms in a double-well confinement potential. At low temperatures (typically ~ 100 nK) the atomic gas forms a BEC 
. _ state, which is a many-particle product state Y\^[AxL{Y^) + BxRi^i)], where Xl/-r(i") is a single particle state localized 
3 . in the left/right well. This wave function is clearly not suitable to define a microscopic two state system — it describes a 
■ " ' collection of non-interacting indistinguishable microscopic (single particle) two state systems that have been discussed 

earlier. A similar system with repulsive scattering favors a Mott state of type 7^[ni=i XL{i^i)\\Wi=N/2+iXR{^i)], 
where V denotes permutation of particle indexes. A schrodinger-cat state ^niXi(''i) + BW^XRi^i) appears only 
in the limit of sufficiently strong attractive interaction. Such a two-state system, however, is difficult to handle. 
The transitions between the left and right macroscopic states are greatly suppressed. The tunneling exponent is 
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proportional to the number of particles. Coherent evolution appears on reasonable time scales only in the limit of 
vanishingly small barrier. 

The outlined difficulties can be overcome in a two-state system defined by macroscopic states that span the same 
spatial region, e.g. supcrfluid (BEC) current states. The simplest geometry featuring such persistent current states 
is a ring. It is inspired by the celebrated superconducting flux qubit |16l llSl Il9| — a superconductor ring with a 
Josephson junction. In a superconductor ring system the qubit is encoded by the magnetic flux states which are due 
to the supercurrent of charged particles (Cooper pairs of charge 2e). Low energy dynamics of these states is due to 
an interplay of three energies: the magnetic inductive energy of the ring, the electric field charging energy, and the 
energy associated with the Josephson junction — a cut in the superconducting ring. Neutral atomic BEC confined in a 
ring featuring weak link (Josephson junction) does not have the first two of these energies. Nevertheless, it turns out 
that such BEC system can still mimic the behavior of superconducting Josephson ring device in certain cases [20| . 

In the following we address the essential physics of the BEC Josephson ring system. We derive an effective low 
energy description (Sec.|lll and show that it is suitable to introduce a macroscopic two-state system, or qubit (Sec.lIIip. 
We discuss possible approaches to single and multi-qubit dynamics (Sec. IIV|) and outline measurement procedures 
(Sec.|Vl. We do not discuss the advantages of BEC- Josephson qubit systems over more traditional qubit designs. The 
goal of this investigation is to introduce an approach to quantum information based on trapped BEC persistent-current 
states. While we give some analysis of a BEC-based qubit and discuss basic quantum operations, we do not perform 
a comprehensive study of single- and multi-qubit gates. Rather, we outline basic principles and state challenges that 
will have to be addressed to fully understand the potentials of this system. 

II. BEC-JOSEPHSON SYSTEM 

In this section we present a detailed review of a single persistent-current BEC-Josephson device. The system is 
inspired by celebrated superconducting Josephson ring — the simplest geometry of a superconductor flux qubit [16| . 
The starting point for our system is a toroidal BEC |2l| - [23| — a BEC of neutral atoms confined in a ring-shaped 
potential. Since the atomic system is in a Bose-condensed state it maintains phase coherence along the ring. As 
the result, the phase fiux through the ring can change only by multiples of 27r, i.e. 27ri^. The integer i^ is a winding 
number. It defines quantized current states: recall that the superfluid current is proportional to the gradient of phase 
along the ring. To have a non-integer winding number the phase has to jump discontinuously at some point along 
the ring (phase slip). This is possible only if the BEC density at that point vanishes. The energy cost of creating 
vanishing density, or a node, comes from kinetic energy and is usually high. However, if a sufficiently high potential 
barrier is introduced, the node can be easily created inside the barrier. As the result, the phase across the barrier can 
change only by a fraction of 27r, causing a small increase in energy, ~ Ej (Josephson energy). The interplay of this 
energy with the energy of the BEC in the bulk of the ring creates rich potential landscape which can be exploited to 
create and manipulate persistent current states. In the following we will derive an effective low-energy action suitable 
to capture essential physics of persistent current states. The derivations will largely follow Ref. [20|. 

The system of interest is shown in Fig. [TJ see also Ref jlj . As soon as atoms are cooled down to nK temperatures 
the trapping is achieved by two laser breams. With the appropriate detuning they can be set effectively attractive 
(potential is lower in region with higher light intensity). The main beam is directed horizontally, see Fig. [TJa), and 
focused to a few /xm in the BEC region. It oscillates to the sides with frequency ~ 100 Hz. The atoms move much slower 
and observe only a time-averaged beam. A virtually fiat 2D-plane potential well is created. The other "painting" 
beam scans the 2D-sheet trap from above creating time-averaged effective attractive or repulsive potential in addition 
to the primary 2D trap. The potential profile corresponding to the BEC-Josephson system is (schematically) shown 
in Fig. [Ijb). It is a toroidal potential well of radius 10 — 100/xm (circumference L ~ 60 — 600/zm) and thickness of 
about 1 — 10/im. The 1/im-size potential barrier is placed at some point of the toroidal well. The barrier can rotate 
along the ring with variable frequency (typically ~ Hz) to stir the condensate. 

The trap with parameters outlined above permits several regimes of BEC dynamics: (a) Complex three-dimensional 
(3D) flow. If the thickness of the BEC ring is significant compared to its healing length, ^ , vortex excitations can 
form in the bulk of the BEC ring. This regime has been explored in [2J]. (b) Quasi-one-dimensional (ID) flow along 
a large-circumference ring. This regime is realized when vortex excitations are squeezed to higher energies by tight 
(of the order of few ^) transverse confinement, while the circumference of the ring is large so that the kinetic energy 
of atonuc motion (BEC flow) corresponding to i' ~ 1 is negligible, (c) Quasi-ID small ring limit. In this case kinetic 
energy of rotation plays an important role. We will focus on the quasi-lD regime of BEC dynamics, i.e. (b) and (c). 

We start with the microscopic description of a quasi-lD BEC system and investigate the partition function, Z, at 
low temperatures. 
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FIG. 1: (a) Essential components of dynamically manipulated "painted potential" trap T]. (b) Ring-shaped BEC (schemat- 
ically) formed in the toroidal potential well cut with the Josephson barrier. The effective (time-averaged) potential profile is 
created by the "painting" beam on top of 2D-sheet confinement. 



where ip is the boson field variable. Our goal is to formulate the effective action in terms of the phase difference across 
the junction (p = ip{Q) — (p{L) integrating out other degrees of freedom. All the parameters are assumed in the units 
of a quasi-lD system; we also set H = 1. The microscopic action S is 



S = So + Sj, 



(2) 
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Sj = 



drj [^*(0, T)^(i, r) + r (L, t)V'(0, t)] . 
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Here P = l/ksT, and the position along the ring, a;, is measured from the junction edge. It is convenient to introduce 
parametrization ipix^r) = ^ p -I- <5p(a;, r)e"''^^''^\ where p is the mean-field BEC density. This density is constant, 
p = N/L, everywhere along the ring except near the Josephson barrier where it approaches zero. Spatial variations 
of density in this region are lumped into the effective Josephson constant J'. With these definitions we obtain the 
standard Josephson action 
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where Ej = JN / L is Josephson energy. Calculations for the bulk action So are more involved. We first rewrite So 
in terms of fluctuation fields 
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The low-energy physics is governed by long- wave fiuctuations I20|,l25[, i.e. V-terms and also (f> are considered small. 
Another small parameter is the magnitude of the density fluctuations. Therefore it is appropriate to keep only the flrst 
three terms under the integral sign in Eq. ^. After this approximation the density fluctuations enter in quadratic 
form, Z ~ j VSpeyL^l—ASp^], and can be integrated out. We obtain 



So = S, 



MF 



is: 
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[c'fc2+w2](p2(fc,a;), 



(7) 



where the velocity of sound is c^ — pXb/tti- This is the standard phonon action formulated in terms of phase 
fluctuations (density and phase fluctuations are conjugate variables). We will focus on the tractable limit for which 
sound fluctuations along the ring are bounded by the Josephson junction and, hence, quantized by the boundary 
condition 92,(^(0, r) = dxip{L,T) = 0. This is appropriate as soon as the level spacing min(c/c) = 27rc/L in the 
excitation energy spectrum is much larger than the Josephson energy per particle Ej /N, i.e. c~^ J . In this limit 




FIG. 2: Typical effective potential energy profiles (in units of Ej) as functions of the phase across the Josephson barrier. 
For comparison, pure Josephson energy contribution is given, (a) Josephson energy, Ej, is small compared to effective the 
"inductive" energy (kinetic energy of rotating BEG), (b) Ej is comparable to the effective "inductive" energy, (c) Same as (a) 
but with rotating Josephson barrier (in rotating frame of reference). 



the action (O can be reduced to a local action by integrating out ip'^{k,u!) assuming that (f>{uj) = ip{0,uj) + ip{L,uj) 
It is accomplished by introducing the functional (5-function according to 



Z = Y[ I V(f>Vipe-^S[(t){uj) - ip{0, uj) + Lp{L, uj)] 



(8) 



and using the identity 5(a) = {2tt)^^ J dAexp (lAa). The fields ip and z can then be integrated out. Note that in doing 
so we expand field ip{x,T) in cos(7rnx/L), thus satisfying the boundary conditions above. After some calculation we 
notice that the partition function takes the form Z — J 2?(/)exp (— iScff ), where 
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In the limit of large ring, L -^ oo (regime (b) as discussed earlier), we recover the Caldeira-Leggett dissipative action 

MM 



S,e ^ ^ J ^\Ljmtj)\^ ~ Ej J dr cos (bir), 



(10) 



which describes dissipative particle dynamics in a periodic potential. 

We are interested in the "coherent" limit when L is finite and a; — )■ (regime (c) as discussed earlier). In this case 
we obtain 
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Note that all the terms in Eq. ([TT|) resemble those of the superconducting flux qubit. In the latter system, however, 
the first two terms originate from classical electromagnetism [16]: the (/)-term, (fi,^C/8e^)0^, is due to the electric field 
of the capacitor formed by different sides of the Josephson junction (and the ring); the (j>^-term., (?i^/8e^£)0^, is the 
result of the magnetic field energy induced by the supercurrent. In our case the kinetic term is the consequence of 
phonon quantization, while the 0^-term is due to the kinetic energy of the collective motion of the particles — superfiuid 
current. 

Finally, we consider the case of a rotating Josephson barrier. At this point it is already clear that the effect of 
rotation should contribute to the 0^-term of the above action. As we will see shortly, the rotation of the trap is 
indeed equivalent to a flux of external magnetic fleld used to tune the superconducting Josephson rings. Consider the 
rotation [13| of the form V{x cos(rit) + y sin fit, x cos(rit) — y sin fit, z). In the rotating frame of reference x ^- x + vt, 
V = Q,L/2iT. This change adds the term J dTpv(j){T) to the above action. Introducing a = vnJL we obtain 
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where 0o = mL^fl/2TT. The parameter a, together with (J)q, determines the shape of the effective energy, i.e. the 
last two terms under the integral sign. Three situations are possible, see Fig. [2 (a) When the tunneling through the 
Josephson barrier is negligible (a — >■ 0) the atomic system condenses into the rotating BEC state defined by (f>o. (b) 
When the tunneling becomes significant the stationary system ($7 = 0) develops a set of metastable current-carrying 
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FIG. 3: Typical double-well configuration: three lowest energies and two basis qubit states are shown. 

states (local minima of the effective potential). The metastable states appear when both the first and the second 
derivative of the effective potential with respect to (j) vanish. The corresponding set of equations is 



tan( 



= — l/cos( 



(13) 



The first pair of metastable states appear when a > ai « 4.6 (0s, i ~ ±1.437r), the second when a > a2 ~ 10.95, and 
so on. (c) Finally, when rotation ^q ~ tt (winding number 1/2) and a > 1 a double-well potential forms, see Fi^[21 
In this review we are interested in the latter case. The analysis of the previous two situations can be found in [20i |. 

III. BEC CURRENT-BASED TWO-SATES SYSTEM: BEC QUBIT 

A. Formulation of a two-state system 

Consider the double- well regime with ipo ^ tt, a > 1. Starting with the classical single-particle action (fT2)) 
analytically continued to real times, we follow standard quantization procedure, see [16]. While such procedure leads 
to the correct result, the reader should refer to a more detailed derivation based on the analysis of the many-body 
wave function 28] . The Lagrangian corresponding to the real-time action is 



^^^^(tf-Ejum)), 
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The motion of the effective "collective" particle can be formarly quantized as follows. First, we introduce canonical 
momentum P and perform a Legendre transformation to formulate a classical Hamiltonian description 



p = dc/d,p^^m, 



'H = P(j)-C. 



We obtain 
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where /i = arj^ /12, rj = ^J p/m\. The quantization is performed introducing P ^ P = —d/d(f>, so that [P,(f>] = 1. 
Finally, we arrive at the effective Schrodinger equation describing the distribution of the collective variable (f> 



Eij{^) = EjHiji^), 



^-i(-^) +^(^)- 



(18) 



The above quantization procedure invokes a conceptual question [19|: The quantization is done on the quantity 
which is defined by purely quantum processes, such as BEC condensation, quantized phonon vibrations, and Josephson 
tunneling. Is it appropriate to use such quantum approach to (J14p at all? A similarproblem with quantization of 
the collective variable takes place in the standard superconducting Josephson ring [16|. In the latter case (due to 
complexity of the microscopic description) this question is resolved primarily by experiment. In our case it can be 



done by constructing an appropriate many-body wave function. It turns out that the superposition of condensates, 
each corresponding to a certain value of phase (f>, weighted by ipi'P)^ reproduces Eqs. ([T51) exactly in the limit of large 



number of particles 1/A^ <C l/rjy/a — 1 ^ VN/ij; further limitations on r] should be enforced to maintain acceptable 
signal-to-noise ratio during measurement [28||. 

The effective potential U and the corresponding energy levels are shown in Fig. [3] When (Jjq '^ n the potential 
acquires double well shape. The splitting energy, 2e, between the two lowest energies is controlled by the height of 
the barrier ~ (a — 1)^/^/q;'^/^ and effective mass fi. For a — 1 ^ 1 and r; ~ 10 — 100 the energy e can be set an order 
of magnitude smaller compared to the energy gap to the third level and a two state approximation can be used. In 
the case oi cfio — tt we obtain 

i\^)=Ho\i^), Ho=ea,, (19) 

where state 1-0) — ^ (0, 1)"^ corresponds to the symmetric ground state of ^TE\\ and 1-0) — ?► (1,0)"'" is the first excited 
(antisymmetric) state. The splitting energy e can be easily calculated numerically. An analytical estimate is e ^ 

(c/L) exp[-2V6r?(a - l)^/^/^]. _^^__ 

Finally, we should point out that the coherent description (J14I18I19P relies on the effective low-energy action 
(|lip . It was obtained assuming uL/c ^1. At sufficiently large w dissipative corrections become important and 
the "coherent" description p^ is no longer valid. In the case of the double- well configuration ((/)o = tt) w can be 
estimated as the frequency of an instanton oscillating inside the barrier (which is potential well for the instanton) . We 
obtain w ^ Ej^/a — l/arj — c^/a — 1/L. Therefore the coherent description is appropriate for small barriers when 
a — 1 ^ 1. On the other hand, a — 1 should not be too small so that different </> states are still distinguishable during 
measurement. 

B. Single-qubit operations and initialization 

Unlike in solid-state quantum computing systems, we cannot use strong measurement procedure [29| to initialize 
the BEC-qubit directly. The BEC-Josephson two state system is based upon interplay of different persistent-current 
states. At the same time, measurement of BEG currents is usually destructive, as we will see in the last section. 
Hence, we should resort to an indirect measurement using a heat bath. Fortunately such measurement is inherent to 
the BEC cooling process — evaporative cooling: this process involves lowering the trap potential to release hot "vapor" 
particles; the remainder of the gas re-thermalizes due to scattering. The preparation starts with cooling of the trapped 
atomic gas. Depending on target winding number, two strategies exist: (i) For winding numbers < 1/2 the pure ring 
potential should be used; no additional (Josephson) barrier is initially required, (ii) For winding numbers > 1/2, a 
high impenetrable barrier should be set rotating with the frequency fJ > fli/2 = tt^/ttiL^. 

We will focus on the case (i). Upon cooling the zero-current (0 ~ 0) condensate forms. Due to the effective 
kinetic term the condensate is in a superposition of different 0-states distributed in the vicinity of ~ with A0 
determined by the effective mass /x. As soon as the Josephson barrier is raised and set rotating by gradually increasing 
the frequency of rotation to satisfy cpo ~ n, this distribution will be adiabatically moved to form the ground state 
of the double-well configuration. Fig. [31 The key is to keep control over the angular acceleration, i.e. change of H. 
with time. In order to have the ground state fully occupied at all times one should keep dU{4))/dt <C (£" — -Eq), or 
00 <C [El — Eq), where Ei — Eq is the gap between the first two energy levels of Eq. (TTS)) far from the double well 
configuration. We obtain ■mL^il/2TT ^ c/L. In the vicinity of the symmetric double- well configuration (0o — tt), the 
Landau-Zener mechanism [30] can be used to drive the system to the antisymmetric state (similar to a single particle 
double- well system). 

To address single qubit operations consider the system close to the symmetric double well configuration with 0o ^ i"- 
In the basis of the symmetric and anti-symmetric states the Hamiltonian takes the form 

H ^ £ay_^ a^, (20) 

a 

where (0)oi is the off-diagonal matrix element of in the same basis. Both phase and flip [23] operation can be 
carried out by tuning the parameters e and 0o (i.e. the Josephson tunneling and frequency of rotation of the barrier). 

IV. MULTI-QUBIT DYNAMICS 

It is well known that superconducting flux qubits can be easily coupled via the magnetic field [18] generated by 
supercurrent in each qubit. However, all the currents in the cold- atom-based qubit are flows of neutral particle. Hence, 




FIG. 4: Stack arrangement; BEC-Josephson rings are placed one on top of the other, so that small tunneling is present between 
the rings. The energy of the system is minimized when the rings are in phase. 

they do not interact with electromagnetic fields. The only interaction present in the system is contact scattering. 
Therefore BEC-Josephson rings must come into contact with one-another or with some common (coherent) BEC 
media. Below we discuss the simplest stack arrangement based on this idea. While more advanced schemes are 
possible, we will focus on demonstrating the basic principles necessary to construct such an interaction in the simplest 
possible model. Multi-plane geometry, however, is not easy to implement with more than a few planes — investigation 
of more complex (in-plane) coupling schemes might be necessary to address the issue of scalability. 

A. Interaction between a few BEC rings 

Consider a stack of M similar BEC-Josephson rings placed one on top of the other in z direction (perpendicular 
to the 2D confinement planes) with the common axis of rotation and the same rotation frequency. The separation 
between the rings is such that weak tunneling is present between the equivalent parts of the adjacent rings, see Fig. ID 
The Euclidian action corresponding to such configuration can be cast in the form 
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The three components are: (i) The effective ID action describing the bulk of each BEC ring 
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(ii) The Josephson tunneling contribution for each ring in a stack 
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(iii) The tunneling between adjacent rings in the stack 
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Here B is the tunneling amplitude between adjacent rings. As before we assume that both Ej = J p and Eb = Bp 
(Josephson energy due to the tunneling between the rings) are small compared to the quantized phonon excitation 
energies, min(cfc) = 27rc/L, in the bulk. We will also assume that all the rings in the stack are arranged such that the 
Josephson barriers are on top of each other. In such system the phases at the same point, x, along the rings do not 
change significantly between different rings 1^9^(2;, r) ~ 1^9^+1(2;, r) — the variation of (pi for each qubit state is relatively 
small {(pi ~ tt) and ipi{x,T) w (pix/L with small fluctuations around this value. We can use the hydrodynamic 
parametrization as before ipi{x,T) = ^ p + Spi{x, T)e'^'^^^^''^' and, to the leading order in fluctuations, obtain 
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Here (and in the following) we suppress the arguments of the fields and the limits of the integration to shorten 
notation. We also expand Sn to the second order in ipi{x, t) — ipi+i{x, r). Higher-order corrections are irrelevant since 
Eb is a small quantity itself. Integrating out the density fluctuations for each ring as before, we obtain 
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Further integration is more involved but follows the same steps as in the case of the single BEC- Josephson ring: (i) the 
phases (/)i across the junction are introduced inserting S[(j)i{T) — {ipi{0, t) — ipi{L, r)}]; (ii) the (5-functions are expanded 
using the identity S{ai) = (27r)^^ / dAi exp {iAitti); (iii) finally, the fields ipi and then A^ are integrated out (Gaussian 
integrals) in the same manner as before. We should note that the Gaussian integrals over ipi are not diagonal in the 
index that counts different rings, hence, a diagonalization has to be done to perform the integration. Fortunately we 
do not have to compute these integrals since they give only a constant shift to the free energy (action) of the system 
which is not important here. As soon as phonon and auxiliary fields are integrated out the second sum of action (j26p 
turns into 
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Following Sec. [IT] we keep only the leading frequency terms. Hence, only the zeroth order in ujL/c should be retained 
in K2{iL'). After some algebra we obtain the effective low-energy action describing the entire stack 



'-'eff — 






2lx'*' + 2£l'*-+''"*' + '^-'™* + iI^''' 



(/),+l)^(l -Si^M) 



(30) 



The form of this action is not surprising. Indeed, by allowing the tunneling between the rings we increase the energy 
of the entire system if those rings are not in phase. This is manifested (to the leading order in the phase difference) 
by the appearance of the last quadratic term in the action, see Eq. (j30p . 



B. Interacting qubits 



We have demonstrated that a set of M BEC- Josephson qubits can interact with the total effective energy as a 
function of the phase-slip phase difference given by 



t^(01,02 



M 

E 



I 

2-J 



bof - cos( 



J 



E( 



(31) 



However, it is still necessary to prove that potential (1311) is sufficient to implement a universal set of gates [29|. Below 
we demonstrate that qubit-qubit coupling given by (I3ip together with a single qubit rotation discussed earlier can be 
used to implement any two-qubit gate. 

Consider a stack of only two BEC- Josephson rings. In the basis of symmetric/antisymmetric states the Hamiltonian 
becomes 



H = Hi+H2 + 



m^)ii 1 



H, 



£(7, 



J 

-K 2tiB 



{4>)oiol 



We tune e — >■ and 



a J 

■n — I'KaB I J . In this case the evolution operator (in real time) is 



U{t) = u'JMt 



exp 



.2tm]n 1 2 






(32) 
(33) 

(34) 



Apart from single qubit rotation Uxa- cTx ^ o' = {<^XT<^yT<^z} this is a natural representation 31] of a (SWAP)^ gate 
with f = 2tB {(j))ffi / Cf . This gate, in combination with single-qubit rotations, is sufficient to implement any two-qubit 
gate [3l|. For instance, one CNOT gate can be realized with the help of two (SWAP)^/2 and single-qubit gates j32l] . 
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FIG. 5: Calculated Time-of-Flight (TOF) measurement outcomes. During the TOF measurement BEC-Josephson qubit wave 
function should collapse to a certain (p state. The density distribution of that state can be captured by absorbtion imaging 
during the TOF expansion. The two-way arrow indicates typical range of (j) expected for measurable two-state BEC-Josephson 
system. 
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FIG. 6: Time-of-Flight density (calculated) at the center as a function of the phase-slip (j) of the ring system before the 
expansion. A typical effective potential corresponding to BEC-Josephson qubit configuration is shown for comparison {(I>l/r 
marks the average position of states in the effective two state system). 



MEASUREMENT 



The only reliable mechanism to measure a cold atom system is via observation of its density. Moreover the atomic 
cloud has to undergo significant expansion before images with suflficient spacial resolution can be obtained. The 
latter is usually accomplished by Time-of-Flight (TOF) experiments [l^l : the trapping potential is turned off and the 
density of the atomic cloud is measured after some time of free expansion. As soon as the confinement potential is 
off, the system is just a collection of atoms that will fly away from the position of the trap, each with the momentum 
supplied by the confinement and motion within the trap. Due to the tight confinement, the released atoms acquire 
substantial velocities which triggers fast expansion and decrease of density. The scattering between particles can be 
neglected, since it usually leads to effects on a much longer time scale. The wave function of the free-expanding BEG 
is 



*(ri,r2,...,i) 



(27r)3 (27r)3 



..*(pi,P2,...)e 



irp — i{p /2m)t 



(35) 



At sufficiently large times only terms with r ^ pt/2m contribute to the integral. As the result, TOF measurement 
provides an image of the density in momentum space. This is ideal to capture the distribution over different current 
states present in (|18p . Consider different moments of density in the momentum space m„ = (/5(qi).../5(q„)), where 
/5(q) = J drdr'e''i'^''-'''^^r')^P{r), 



m„ 



Vi>*V'4} / dridr;...dr„dr;e''ii(""i-'"'i)...e*'i"(''"~''"V*(r'i,0)V'(ri,0)...V'*(r;,0)i/;(r„,0)e- 



(36) 



Here S is given by Eqs. ([3|) and ([4]). The expansion in terms of small fluctuations ^(r, r) — x</)(T)(r) + (5i/'(r, r) yields 

dcpP{4>)\xMi)?-\xMn)? + 5'[/(qi, q2) + /(q2,q3) + ...] + 0(.9'), (37) 



where the (path) integrals over (j) with r 7^ has been lumped into P{<j)). The function /(q, q') is a combination of two 
particle Green's functions describing scattering processes in and out of the condensate, see [33j. These two-particle 
correlations are suppressed il3] by the gas parameter r^ g, which limits the resolution of the measurement. For the 
range of parameters appropriate for the macroscopic two state BEG-Josephson system the gas parameter is small. 
Apart from the g^ correction, moments (|38)) define a stochastic process: P{(f>) plays the role of the probability to find 
the system at state cj) as the result of the measurements. Hence, by measuring the density distribution in momentum 
space (via TOF) we measure the value of cf). At each experimental run a different density distribution is expected 
with the probability determined by jV'l^)!^. 
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The possible TOF measurement outcomes are given in Fig [5l where we plot the momentum-space density dis- 
tributions |x0(q)P for several values of from to 27r. Different (p states are clearly distinguishable. For a more 
quantitative judgment we propose to measure the TOF density at the center of the trap. The central TOF density, 
|X0(O)P = l/'^rx0(r)p, depends on the phase i^(r) in a straightforward way. In the bulk of the ring this phase is 
linear (in this case kinetic energy is minimized) (p{r) = <j)x/L and, hence, X0(r) ~ g^<t>^/L_ ^g ^-j^g result we obtain 

fn\ 'l'=o ^„^ sin^(0/2) 
Ptof{0) = PtofW ,ii2\2 ■ (3^) 

For the variation of (j) relevant to the BEC-Josephson qubit system, see Fig [SI this is an approximately linear function 
4/7r2-8(0-7r)/7r2. 

In this review we have considered only the very basic properties of the cold atom BEC-Josephson rings relevant 
to quantum information. A detailed investigation is necessary to assess noise effects as well as scalability. A simple 
proof-of-principle experiment with a single BEC-Josephson macroscopic two-state system should also be helpful in 
determining the direction for further theoretical development. 
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